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Abstract 

In this paper, an evolutionary strategy (ES) method is introduced as an optimization 
approach to solve problems in the manufacturing area. The ES method is applied 
to a case study for milling operations. The results show that it can effectively yield 
good results. 
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Insert Nomenclature here 


1 Introduction 


In manufacturing industry, it is important to determine the optimal machining 
parameters in order to maximize total proht rate and to increase the quality 
of the hnal product for machining operations. Traditionally, the selection of 
parameters is carried out by the experience of planners or with the help of 
machining data handbooks. These may not guarantee the optimum perfor¬ 
mance and minimization of costs. Therefore, a number of researchers have 
tried to deal with the optimization of machining parameters using different 
approaches. Compared with some deterministic methods [1-4] evolutionary al¬ 
gorithms (EAs) are more attractive to engineers since EAs are robust, effective 
and easy to implement. 

Many different types of EAs have been proposed for optimizing machining 
parameters of milling operations, such as genetic algorithm (GA) [5], simulated 
annealing (SA) [6], ant colony algorithm [6], immune algorithm [7], and some 
hybrid algorithms [8,24-26]. 

Evolution Strategies(ESs), originally developed by Rechenberg and Schwefel 
[12,13], are algorithms which imitate the principle of natural evolution as 
a method to solve parameter optimization problems [14,27]. ESs can reject 
infeasible individuals directly (also called “death penalty”). It is probably the 
easiest way to handle constraints. Since in the model introduced in sequel, 
the feasible search space constitutes a reasonably large portion of the whole 
search space. This strategy seems feasible. Moreover, it is also computationally 
efficient, because when a certain solution violates a constraint, it is assigned a 
htness of zero. Therefore, no further calculations are necessary to estimate the 
degree of infeasibility of such a solution. To the best knowledge of the authors, 
no one has used ESs to optimize machining parameters of milling operations. 
So in this paper, we use an ES to solve this problem. The results show that it 
can effectively yield good results. 


2 Mathematical model 


Depth of cut, feed rate and cutting speed have the greatest effect on the 
success of the machining operation. Depth of cut is usually predetermined 
by the workpiece geometry and operation sequence. It is recommended to 
machine the workpiece with the required depth in one pass to keep machining 
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time and cost low, when possible. Therefore the problem of determining the 
machining parameter is reduced to select the proper cutting speed and feed 
rate combination. The mathematical model developed by M. Tolouei-Rad et 
al. [16] is considered in this work. The model is also considered in several 
papers ( [6,8,17,18]). 


2.1 Objective function 


The main focus of this work is to maximize the total proht rate and can be 
determined by 
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The unit time to produce a part in the case of multi-tool milling can be defined 
by 


T, = + y Vr' + E (3) 
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2.2 Constraints 


In practice, possible range of cutting speed and feed rate are limited by the 
following constraints 

1. Maximum machine power 

2. Surface hnish requirement 

3. Maximum cutting force permitted by the rigidity of the tool 

4. Available feed rate and spindle speed on the machine tool 

2.2.1 Power 

The machining parameters should be selected such that maximum machine 
power is used. The required machining power should not exceed available 
motor power. Therefore the power constraint can be written as 

C,Vf-^ < 1, (4) 
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where 
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2.2.2 Surface finish 


The required surface hnish Ra, must not exceed the maximum attainable 
surface hnish Ra{at) under the conditions. Therefore the surface hnish for end 
milling becomes 



Cd < 1, 

(6) 

where 

318[tan(/a) + cot(ca)]“^ 

Ce - P 

■^a(at) 

(7) 

And for end milling 

Cif‘ < 1, 

(8) 

where 

^ 318(4<i)-‘ 

^a{at) 

(9) 

2.2.3 Cutting force 



The total cutting force F/ resulting from the machining operation must not 
exceed the permitted cutting force Fc(per) that the tool can withstand. The 
permitted cutting force for each tool has been considered as its maximum limit 

for cutting forces. Therefore the cutting force constraints becomes 



CsFc < 1, 

(10) 

where 

Cs = 1/Fc(per). 

(11) 


2 . 2.4 Speed limits 

1. Face milling: 60-120 m/min 

2. Corner milling: 40-70m/min 

3. Pocket milling: 40-70 m/min 

4. Slot millingl: 30-50 m/min 

5. Slot milling2: 30-50 m/min 


2.2.5 Feed rate limits 

1. Face milling: 0.05-0.4mm/tooth 

2. Corner milling: 0.05-0.5mm/tooth 
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3. Pocket milling: 0.05-0.5mm/tooth 

4. Slot millingl: 0.05-0.5mm/tooth 

5. Slot milling 2 : 0.05-0.5mm/tooth 


3 The ES method 

3.1 Introduction of ESs 


Evolution Strategies can be understood as ‘intelligent’ probabilistic search 
algorithms which are based on the evolutionary process of biological organisms 
in nature. 

The procedure of one type of ES can be described as follows. For the fi initial 
individuals, in optimization terms, each individual in the population is en¬ 
coded by a real number which represents a possible solution to a given prob¬ 
lem. Then a recombination procedure and mutation procedure are executed to 
produce new ‘offspring’(i.e. children) solutions with 77(77 > jj) individuals. The 
htness of each individual in ‘offspring’ solutions is evaluated with respect to a 
given objective function. After evaluation procedure the offspring solutions are 
sorted and the last rj—jj, individuals are deleted. This reproduction-evaluation- 
selection cycle is repeated until a satisfactory solution is found. 

A more comprehensive overview of ESs can be found, e.g. in [14,19] and 
references therein. The applications of ES could be found in [21-27] 

The basic steps of the procedure can be shown as: 

Generate an initial population; 
repeat: 

Recombinant and mutate individuals to produce children; 
Evaluate fitness of the children; 

Select the population from the children; 
until a satisfactory solution has been found. 

By means of the above ES procedure, a computing method for optimizing 
machining parameters of milling operations is described in the following parts 
of the section. 
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3.2 Representation and fitness function 


We denote X as a vector composed of two arrays representing cutting speed 
and feed rate, i.e., Xi = V and X2 = f, X = (a;i,a; 2 ). The target function (1) 
is taken as the htness function in the ES. 


3.3 Initial population, recombination and mutation 
3.3.1 Initial population and Recombination 

The initial population consists of/i individuals ^(0) = • • • ,X(^^(0)}, 

where X«(0) = {X\ afi. The initial components of X are generated randomly 
from feasible domains. All initial Uj are valued by 3.0 here. 

Select two individuals: 

(X^cr^) = ((xj,--- ,x])fia\,--- ,a])) 
and (X^ 0 - 2 ) = {{xj, ■■■ ,x‘f), ■■■ , af)). 

There are two types of recombination operators: 

• discrete, where the new offspring is 

with qi and g' equal either 1 or 2; 

• intermediate, where the new offspring is 

(X, cr) = {{axl+{l—a)xl, • • • , {axl+{l—a)x‘f), {aal+{l—a)(7l, • • • , aal+{l—a)(7f)), 
with a G (0,1). 

By the suggestion of Schwefel [13], the discrete recombination operator is 
executed on X, and the intermediate recombination operator is executed on 

cr. 


3.3.2 Mutation 

Apply mutation to the offspring (X, cr), the resulting new offspring (X', cr') is 
obtained, where 


a[ = ai- exp(r' • X(0,1) + r • Xi(0,1)) 
and x'i = Xi + o'i ■ Xj(0,1), 


( 12 ) 
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in which, 

7V(0,1) -a random number reconciled standard normal distribution; 

7Vj(0,1)-a random number reconciled standard normal distribution aimed 

at ith component; 

t' -global coefficient; 

r-local coefficient; 

I -the number of components in X. 

According to the suggestion of Schwefel [13], r' is valued by (\/^)~^ and r is 
valued by (\/2\//)“^ in the computation. 

Notice that X has box constrains. So the offspring generated by mutation 
procedures may not be feasible. Then we just let it take the bound value, i.e. 
if Xi > Xi (or < Xj), let Xi = Xi (or = Xj). 


The recombination and mutation steps will not stop until rj offsprings are 
generated. 


3-4 Evaluation, Selection and Stop criterion 


We evaluate every {X, cr) by its htness function and sort them. If an individual 
is infeasible, the htness value is 0. Then we choose hrst fi individuals as new 
parents. In our ES method, we let p, = 15, rj = 105 and set a variable to record 
the current best htness value. If a htness value is better than the record, we 
will update the record. The algorithm will stop if the record keeps unchanged 
after 1000 iterative loops. 


4 Case Study 


The component as shown in Figure 1 is to be produced using a CNC milling 
machine. It is desirable to hnd the optimum machining parameters, which 
result in the maximum proht rate. Specihcations of the machine, material and 
values of constants are given below. Also, the geometric information on the 
required operations and tools is presented in Tables 1 and 2. 


Insert Figure 1 here 


Constants: 

Sp = $25 
cmat = $0.50 
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Co = $1.45 per min 
Q = $0.45 per min 
tg = 2 min 
tie = 0.5min 

C = 33.98 for HSS tools 
C = 100.05 for carbide tools 
w = 0.28 
Kp = 2.24 
W = 1.1 

n = 0.15 for HSS tools 
n = 0.3 for carbide tool 
g = 0.14 

Machine tool data: Type: Vertical CNC milling machine = 8.5kW, 

e = 95% 

Material data: Quality: 10L50 leaded steel. Hardness = 225 BHN 


Insert Table 1 and Table 2 here. 


We use our ES to solve this case and compare it to some other methods in the 
literature. The results listed in Table 3 show that ES can obtain good results 
similar to hybrid methods. 


Insert Table 3 here. 


5 Conclusion 


In this paper, we used an ES method to optimize machining parameters of 
milling operations. The results showed that it can effectively give good results 
and it can be a good alternative in similar problems in engineering. 
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Nomenclature 


O'rad 


ca 

Ci ii = l,...,8) 


Q? Co 

Cm? ^mat-) Ct 

Cu 

d 


e 

F 


f 

Fc, Fc(per) 


G,g 


K 

Kiii = 1,2,3) 

Kp 

la 

m 

N 


n 


P, Pm 

Pr 

R 

Pai Pa(at) 

Sp 

T, Tu 


tmt ts7 Pc 


V, ^hb, Vopt 


W 


W 


z 


Axial depth of cut, radial depth of cut (mm) 

Clearance angle of the tool (degrees) 

Coefficients carrying constants values 
Labour cost, overhead cost ($/min) 

Machining cost, cost of raw material per part, cost of a cutting tool ($) 

Unit cost ($) 

Cutter diameter (mm) 

Machine tool efficiency factor 
Feed rate (mm/min) 

Feed rate, (mm/tooth) 

Cutting force. Permitted cutting force (N) 

Slenderness ratio, exponent of slenderness ratio. 

Distance to be travelled by the tool to perform the operation (mm) 
Coefficients carrying constant values 
Power constant depending on the workpiece material 
Lead (corner) angle of the tool (degree) 

Number of machining operations required to produce the product 
Spindle speed (rev/min) 

Tool life exponent 

Required power for the operation, motor power (kW) 

Total profit rate ($/min) 

Sale price of the product excluding material, setup and tool changing costs ($) 
Arithmetic value of surface finish, and attainable surface finish (^m) 

Sale price of the product ($) 

Tool life (min). Unit time (min) 

Machining time, set-up time, tool changing time (min) 

Cutting speed, recommended by handbook, optimum (m/min) 

Exponent of chip cross-sectional area 
Tool wear factor 

Number of cutting teeth of the tool 
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Table 1 

Required machining operation 


Operation Number 

Operation type 

Tool number 

a (mm) 

iF(mm) 

Ra(^m) 

1 

Face milling 

1 

10 

450 

2 

2 

Corner milling 

2 

5 

90 

6 

3 

Pocket milling 

2 

10 

450 

5 

4 

Slot milling 

3 

10 

32 

- 

5 

Slot milling 

3 

5 

84 

1 


Table 2 
Tools data 


Tool Number 

Tool type 

Quality 

(i(mm) 


Price($) 

la 

ca 

1 

Face mill 

Carbide 

50 

6 

49.50 

45 

5 

2 

End mill 

HSS 

10 

4 

7.55 

0 

5 

3 

End mill 

HSS 

12 

4 

7.55 

0 

5 


Table 3 

Comparison of the results for milling operation 


Method 

Ctj-Unit cost 

T^j-Unit time 

P^-Profit Rate 

Handbook [20] 

$18.36 

9.40 min 

0.71/min 

Method of feasible direction [16] 

$11.35 

5.48 min 

2.49/min 

Genetic algorithm [17] 

$11.11 

5.22 min 

2.65/min 

Ant colony algorithm [6] 

$10.20 

5.43 min 

2.72/min 

Hybrid particle swarm [18] 

$10.90 

5.05 min 

2.79/min 

Immune algorithm 

$11.08 

5.07 min 

2.75/min 

Hybrid immune algorithm [17] 

$10.91 

5.07 min 

2.79/min 

Hybrid differential evolution algorithm [8] 

$10.90 

5.00 min 

2.82/min 

Evolutionary strategy 

$10.91 

5.00 min 

2.82/min 
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